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Abstract We evaluate the deuteron charge, quadrupole, and magnetic form factors using wave functions 
obtained from chiral effective theory (xET) when the potential includes one-pion exchange, chiral two-pion 
exchange, and genuine contact interactions. We study the manner in which the results for form factors 
behave as the regulator is removed from the xET calculation, and compare co-ordinate- and momentum- 
space approaches. We show that, for both the LO and NNLO chiral potential, results obtained by imposing 
boundary conditions in co-ordinate space at r = are equivalent to the /l — > oo limit of momentum-space 
calculations. The regulator-independent predictions for deuteron form factors that result from taking the 
yl — > oo limit using the LO xET potential are in reasonable agreement with data up to momentum transfers 
of order 600 MeV, provided that phenomenological information for nucleon structure is employed. In this 
^ ' range the use of the NNLO xET potential results in only small changes to the LO predictions, and it 

improves the description of the zero of the charge form factor. 

00 

^ : PACS. 12.39.Fe Chiral Lagrangians - 25.30.Bf Elastic electron scattering - 21.45.+V Few-body systems 

^ " 1 Introduction 

O . 

• • Elastic electron scattering from the deuterium nucleus, the deuteron, has long been used as a means to assess the 
. ^ reliability of different nuclear forces. For recent reviews of experimental and theoretical work on the elastic electron- 
: deuteron reaction see Refs. [TJ[2l[3]. In this reaction the ability to vary the momentum transfer, to the nucleus 
^ ■ means the electromagnetic structure of the deuteron can be probed on a variety of different scales. Indeed, in the 
' "impulse approximation" — where two- nucleon currents are neglected — the matter distribution inside the nucleus is 
imaged, thereby providing a direct test of the NN interaction used to predict it. In this paper we examine elastic 
electron-deuteron scattering for q < 1 GeV, where we anticipate that the deuteron's electromagnetic structure will be 
governed by the spontaneously and explicitly broken chiral symmetry of QCD. 

This feature of QCD plays a role even though quark degrees of freedom are not resolved at this momentum scale. 
Indeed, confinement guarantees that for g < 1 GeV it is useful to employ an effective field theory (EFT) in which 
nucleons and pions (and possibly baryon resonances) are the explicit degrees of freedom. Here we will use chiral 
perturbation theory (^PT) which is an EFT in which the nucleons and pions interact with each other in a manner 
consistent with the chiral symmetry of QCD and the pattern of its breaking. xPT does not aim to describe hadronic and 
nuclear structure down to arbitrarily short distances, but is based on the inability of low-momentum-transfer probes 
to resolve the details of short-range interactions. It is thus ideally suited to the treatment of deuteron electromagnetic 
structure in the regime q <1 GeV. 

Chiral perturbation theory organizes the Lagrangian that describes nucleon and pion (and photon) interactions in 
powers of the small expansion parameter P = ^"^^ , where A^sb ~ ^pt^t^I-k represents the scale of chiral-symmetry 
breaking, which is indicated by the p mass and the pion decay constant, while p is the momentum of the particles 
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involved, and TOtt the pion mass. This theory has had considerable success in describing pion-nucleon and meson- 
meson interactions I4,5j. Weinberg proposed to extend xPT to few-nucleon systems by making a xPT expansion for 
the nucleon-nucleon interaction, V , and then taking the non-perturbative character of nuclear systems into account by 
solving the Schrodinger equation based on this potential [6l[7l[8]. The result is quantum mechanics with fixed particle 
number for few-nucleon systems in which the potential (and other operators) has a expansion. This has been 
dubbed chiral effective theory (xET). This approach is reviewed in Refs. [9[ fT0l[Tn[T2] . 

Recently there has been considerable controversy about the correct way to implement this approach [T3lfT4lfT5l 
I16j . It has been argued that Weinberg's original power counting needs to be modified for the short-range or contact 
interaction piece of V if one wishes to obtain regulator-independent predictions over a wide range of cutoffs, A. 
However, as shown already in Ref. pT], in the '^Si-'^Di channel the leading-order chiral potential can be renormalized 
using the single contact interaction it contains. The limit yl ^ oo can be taken, as has since been confirmed by several 
authors p!4l[T8l[T6lll9j . Thus the discussions of Ref. [T3lll4l[T5lfT6] are not relevant to our results for leading-order (LO) 
xET, and we will not contribute to that discussion here. 

However, the situation is not as clear once sub-leading corrections are included in the NN interaction V . In 
particular, while one naively expects that these corrections to V can be treated in perturbation theory, previous works 
have included them non-perturbatively, by solving the Schrodinger equation with the full V dX a, fixed order in xPT. 
In such a calculation the sub-leading contributions to V dominate over the leading-order part of the xPT potential for 
r <C l/m-TT, and the singularity of V at short distances therefore becomes more severe as the chiral order to which V 
is computed is increased. In what follows we adopt such a non-perturbative treatment of V , even though it leads to a 
paradox when the next-to-leading-order [0(P^)] V is employed. In standard xPT, where only nucleons and pions are 
taken as the relevant degrees of freedom, the next-to-leading order (NLO) V calculated for r ~ 1/m^ is repulsive for 
r <C As is well known from the theory of singular potentials [20], this means that the spectrum of the potential 

is predicted in the limit that the cutoff on the potential is removed, i.e. that the xPT potential is taken to be valid 
for all r. In consequence in the limit yl — > oo there is no opportunity to tune the contact interaction and thereby 
reproduce the binding energy of the shallow NN bound state in the '^Si-^Di channel [2T] . 

Therefore, we can only present results for the LO [0(P°)] and next-to-next-to-leading order [NNLO=0(P^)] 
chiral potentials V . Both of those potentials are attractive for r <^ l/m^, which means that the inclusion of contact 
interaction(s) is mandatory if cutoff-independent predictions are to be obtained. Moreover, at these orders, and in the 
^Si-^Di channel of NN scattering, there is no problem with Weinberg's original proposal: the three short-distance 
A^A^ operators his xPT expansion of V predicts should be present at NNLO are sufficient to renormalize the attractive 
long-range potential generated by pion-nucleon dynamics at r ~ l/m^r. 

In Section [2] we compute the deuteron wave functions obtained from the LO and NNLO chiral potentials V. We 
show that when we compute using a momentum-space cutoff on V but demand that the deuteron binding energy be 
reproduced we obtain cutoff-independent predictions for the radial wave functions u{r) and w(r) for r ^ In the 
limit yl — > oo the result of such a momentum-space calculation then agrees with that obtained via the co-ordinate- 
space approach advocated in Refs. [T51l21| . We use existing wave functions to make this demonstration for the LO 
xPT potential V, and we also obtain, for the first time, NNLO momentum-space wave functions in xET for yl > 1 
GeV. This allows us to make the connection between momentum- and co-ordinate-space not just at LO, but also at 
NNLO, and, by extension, for any order at which V is attractive at short distances. 

These wave functions, while not directly observable, do enter in the matrix elements of nuclear current operators 
that are probed in electron scattering from the deuterium nucleus. In particular, to the extent that the impulse 
approximation is valid, the form factors that determine all elastic electron-deuteron scattering observables in the one- 
photon-exchange approximation can be expressed as Bessel transforms of bilinears of these wave functions. We present 
results for these form factors in Sec. [3l These results are based on a chiral expansion of the deuteron current operators 
mug. Up to relative order there are no two-nucleon contributions in either the charge operator, J'^ or the vector- 
part of the current, J, and the impulse approximation holds. The deuteron form factors have already been calculated 
in such an approach (including two-body currents) in Refs. p4 l f25 l [26 j l27j. Here, we follow 26J and focus on tests of 

deuteron structure by computing the ratio G'c/G^'', where the denominator is the nucleon isoscalar form factor. We 
show that cutoff artifacts disappear in this ratio as yl — > oo, allowing us to obtain the first regulator-independent, 

(s) 

LO, results for Gc in the literature. We find that if we use phenomenological input for these results are in good 
agreement with data for Ge- 
lt xET is to be a systematic expansion for A^A^ interactions we should be able to also obtain results at higher 
orders in the P expansion. Therefore in Sec. [4] we compute the results found for form factors when Gc, Gq, and Gm 
are computed with deuteron wave functions found using the NNLO xPT potential. This potential includes various 
contributions from two-pion-exchange mechanisms. Once again we find that the cutoff artifacts in observables vanish 
as yl ^ GO. Importantly, we find that the shift from LO to NNLO results for the form factors is small for q < 600 
MeV. 

This suggests that if we could do a NLO calculation of deuteron structure it would be close to the regulator- 
independent LO result in this kinematic domain. In the absence of a calculation at NLO with yl — > oo in the theory 
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with only explicit nucleoli and pion degrees of freedom this is purely supposition. Nevertheless, we believe that our 
results with the NNLO V and yl — > oo give realistic estimates for two-pion-exchange contributions to the form factors. 
We hold this view because the 0{P^) pieces of V that make the difference between an NNLO and NLO calculation 
are numerically much more important than NLO pieces [28,29,30^. This is in large part because our NNLO V includes 
the most important contributions of the Z\(1232) to the NN interaction p51l5T[|5n] . In particular, if a consistent power 
counting for nuclear interactions requires taking the zi(1232) into account as an explicit degree of freedom [28l[3Tll32| 
our NNLO calculation will produce results similar to an NLO calculation in the theory with explicit Deltas. 

The calculation presented in Sec. 2] shows significant sensitivity in the position of Gc's zero to the details of the 
NNLO xPT potential. However, any attempt to use data on the go for which Gcilo) = to pin down the two-pion- 
exchange part of the NN force is clouded by the presence of two-nucleon pieces of Jq at 0{P^) — the same order 
relative to leading as that at which the two-pion-exchange pieces of V first occur. Similar ambiguities bedevil attempts 
to use other features of deuteron electromagnetic form factors to identify aspects of the NN force beyond one-pion 

exchange. We provide an estimate of the impact on our results of O (^^pjj^ and O (^^^ i pieces of Jq in Sec. [H But 
a full calculation of their contribution when yl ^ oo is beyond the scope of this paper. In Sec. [51 we summarize our 
findings and conclude. 



2 xET wave functions: equivalence of co-ordinate- and momentum-space solutions 

xET for few-nucleon systems results when xPT is applied to derive potentials and current operators which can then 
be used in a non-relativistic quantum-mechanical (fixed-particle-number) framework. Corrections to the static-nucleon 
picture may be treated in such a framework, e.g. here we incorporate the 1/M corrections to V as given by [53]. 1/M^ 
pieces of V and the effects of higher-Fock-space (e.g. ttNN) states can also be incorporated, but will not concern us 
significantly here, although the former will be briefly discussed in Sec. [D 



2.1 Regulator-independent results for deuteron wave functions with the leading-order xPT potential 

As first discussed by Weinberg j6l[7] the leading-order potential in the xET reproduces the time-honoured one-pion 
exchange (OPE) potential 



F(°)(p,p') 



9a ^1 ■ (P' ~P) (^2- (p' - p) 



4/2 



(P' - P) 



Cs + Ct CTl ■ (72 



(1) 



We can now solve the Schrodinger equation in its momentum-space form, i.e. the homogeneous Lippmann-Schwinger 
equation: 

d^p' 



LO/A 



Go{p) 



(2) 



where A is the scale at which the potential V is regulated, and G'o(p) — {—Bd — /M)^^ is the (free, center-of-mass 
frame) two-nucleon propagator, with Bd and M denoting the deuteron binding energy and nucleon mass, respectively. 
The OPE contribution is determined through the pion mass m^r — 138 MeV, the axial coupling constant (for which 
we adopt gA = 1.29) and the pion-decay constant — 92.4 MeV. In the triplet channel, only one linear combination 
of the contact interaction parameters, Cs + Ct, is relevant and is adjusted for each cutoff A to reproduce the deuteron 
binding energy B = 2.225 MeV as previously done, e.g., in [33j . We perform the regularization using exponential cutoff 
functions 

8 \ / '8 \ 

y(°)(p,p')exp ' ^ 



Vl°^(p,p') = exp 



"EL 
A^ 



A^ 



(3) 



Alternatively, the interaction can be Fourier-transformed to co-ordinate-space resulting in 



V{v) = 



Cs + CtCi • 172 



(4) 



The Fourier transformation uniquely determines the long-range part of the interaction. (Note, however, that it involves 
a redefinition of the contact interactions to absorb a piece of that is a constant in momentum space.) The deuteron 
wave function in the pn CM system can be represented as 



47rr 



u{r)ap ■ an + {3ap 



• an ■ f - (Tp ■ an) 



A-pn 



(5) 
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with the total spin s = 1 and = 0, ±1 and CTj, and (T„ the Pauh matrices for the proton and the neutron, respectively. 
The functions u{r) and voir) are the radial ^Si and ^Di components of the relative wave function, respectively. 

Even if we consider only r > and ignore the presence of the three-dimensional delta functions, this is a singular 
quantum-mechanical potential. In the absence of any short-distance regulator, the resulting Hamiltonian is unbounded 
from below and the spectrum contains an infinite number of bound states. This happens for any value of the parameters. 
On the other hand, it is known that physically there is only one bound NN state: the deuteron. The presence of these 
unphysical, spurious bound states should not have any impact on low-energy observables, including the form factors 
at low momentum transfers. We will come back to this issue in Section [3] 

To derive the wave function from Eq. ([4]), we use the representation of Eq. ([5]) and find the usual coupled one- 
dimensional differential equations for the radial wave functions u and w in the presence of a tensor potential: 



u"{r) + Us{r)u{r) + Usd{r)w{r) = --j^u{r), 



'w"{r) + Usd{r)u{r) + 



Ua{r) + ^ 



w{r) — — 7^ii;(r), 



The coupled-channel reduced potential [U ~ 2^pnV with 2/ip„ 

U, = Uc, Usd = 2\/2C/t, 

with 



2MpMn/{Mp + Mn)) is given by 
Ud = Uc- 2Ut, 



(6) 



(7) 



Uc = 
Ut = 



16^/2 



r 



16^/2 



1 



rriTrr 



(8) 



provided that r > 0. 

For our co-ordinate-space calculations, the equations ([6]) are solved subject to the boundary conditions 



u{r) 
w{r) 



A 



rjAse' 



1 + A 

7r 



(9) 



for r oo. Here, 7 = \/ MBd is the deuteron wave number. As is a normalization constant and can be chosen to 
guarantee that 

/•oo 

dr(u2(r) + w'^{r)) = 1 , 



(10) 



and 7] is the asymptotic D/S ratio. 

The form of the singularity of the OPE potential at short-distances implies that, for sufficiently small r, the 
components u and w are given by [M1I181I35] : 



Ushort{r) = As- 



Wshort{r) 



As- 













)"[ 







-C2«e-4^\^ + 23/2|C2^|cos Jf^ 



V2C: 



2Re 



+ 2|C2a|cos 4 




(11) 



CiA and CiB, are normalization constants which have been determined in [18 , i? is a new length scale that enters the 
non-perturbative problem, given by i? = If^. When Eqs. dH)-® are solved in Ref. [TS] the phase is determined 
by the boundary condition at r = 0, and so </) is regulator independent, and is a function only of the scales m^, 7, and 
R. 

The wave functions u and w can now be calculated by applying standard techniques to Eqs. ([6])-([8]) and employing 
the long-distance boundary conditions ([9]) [ISj . The asymptotic D-to-S ratio 77 is then determined by the requirement 
that the numerical solution matches the short-distance behavior Eq. (jlip . In practice, this is done by imposing the 
additional boundary condition uir,^ — \f2w{r^ = 0, which is a direct consequence of Eq. (jlip . at a value of R. 
For this purpose we choose = 0.1-0.2 fm [Hj, and find ?yoPE(7) ~ 0.026333. 

The necessity to choose a finite leads to a small numerical uncertainty of the constants C2A^ C2B. quoted in 
Ref. |18] . These uncertainties can, however, be made arbitrarily small. Furthermore, any deuteron matrix element that 
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is finite when computed with the wave functions u and w obtained in this way is insensitive to as long as Tc is taken 
small enough [36]. Therefore, for all practical purposes, the wave functions found by the technique of Ref. [H] are the 
deuteron wave functions found with a strict contact interaction. The contact parameter Cs + Ct has been converted 
to a boundary condition that is imposed at r = [37j . 

For non-singular interactions, it is clear that any solution obtained in momentum space will be the Fourier-transform 
of the corresponding co-ordinate-space solution. The question arises whether this holds for the solution of the singular 
OPE interaction based on Eq. ([2|). In momentum space the solution is only ever determined for a finite cutoff A, and 
it is not a priori clear that the limit yl ^ oo of the momentum-space wave function will, upon taking the Fourier 
transform, lead to the co-ordinate-space wave function found using Eqs. ([6|)"([8|). 

In Fig.[Tl we show that, for any r > 0, the momentum-space solution does indeed approach this co-ordinate-space 
solution for A — + oo. It is reassuring that one obtains equivalent results in both schemes. This confirms that fitting 
the contact interactions in a momentum-space calculation corresponds to imposing proper co-ordinate-space boundary 
conditions and matching to the short-distance form of the wave functions (fTT|) . 

However, the relationship between A and the numerically necessary co-ordinate-space cutoff Tc remains somewhat 
opaque. For theories in which the long-range-part of the potential is absent one can establish that a sharp cut-off A 
and the matching radius Vc for the Schrodinger equation with a boundary condition are related by rc = ^ [38]. But 
here, given our choice of regulator and the presence of a long-range potential (one-pion exchange) this relation is no 
longer valid. Indeed, especially for w{r), the wave function is only converged for distances 

.»^. (12) 

Therefore we conclude that simple arguments for "reasonable" choices of the momentum-space cutoff A at which 
results converge to the yl — > oo limit can be misleading. 



2.2 Regulator-independent results for deuteron wave functions with the NIMLO potential 

In Sec. [3] we use the (co-ordinate- and momentum-space) wave functions we have obtained thus far to predict the 
electromagnetic structure of the deuterium nucleus. However, this represents only a LO calculation in xET, so before 
performing such a calculation of deuterium electromagnetic form factors we will explain how the two-pion exchange 
(TPE) contribution to the NN potential V is incorporated in xET. 

As outlined in the introduction, xPT can be used to expand the long-range (r ~ l/m.,^) part of the NN potential. 
Here we do this up to NNLO. The explicit expressions in momentum- and co-ordinate-space were given in Ref. |29j . 
Note, however, that the momentum-space expressions of Ref. |29| are opposite in sign to ours, since different conventions 
are used. Also, it should be noted that we have used gA = 1-26 for the TPE part of the NNLO potential. 

For the calculations in co-ordinate-space, we follow the formalism introduced by two of us in Ref. 21j where two- 
pion-exchange effects were included non-perturbatively. For the NNLO interaction the short-range potential behaves 
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Table 1. 



("21 

The four different sets of values for the chiral coefficients from C\l^ that we consider in this work. 



Set 


Source 


ci(GeV-') 


C3(GeV-i) 


C4(GeV-') 


Set I 


■kN [42J 


-0.81 


-4.69 


3.40 


Set II 


NN [M, 


-0.76 


-5.08 


4.70 


Set III 


NN [43] 


-0.81 


-3.40 


3.40 


Set IV 


NN [H] 


-0.81 


-3.20 


5.40 



like [lilMlIlD] 



U. 



TPE 



u. 



TPE 



^6 

Rt, 



Ri. 



(13) 



where 



{Rsf 
{Rsdf 



{Rdf 



35i 



(4- 



128/>2 
128/47r2 



3\/2.9i 



(-4 + 



32/47r- 



■{--1 + 



25a 



24c3 - 8c4), 

I6C4), 
2c3-2c4); 



(14) 



and Ci = Mci are the low energy chiral coupUngs (LECs) appearing in ttN scattering [41]. We will present results for 
the sets of Cj parameters shown in Table [TJ which are the same four sets as were used in Ref. [T^. In contrast to OPE, 
this interaction yields two attractive eigenchannels upon diagonalization. Therefore, the short-distance wave functions 
in the ^Si and '^Di channels take the form (21j: 



Ushort{r) = {C+ 



r 



cos U sm 



Wshort (r) 



2 r2 ^ 

i Rl 
272" 



C- 



r 



sm u sm 







^-^^ COS sin 







1^ 

2 r2 ^ 

1 Rl 
2~ 



(15) 



where the coefficients C+ and C_ and the phases 0+ and (p- are arbitrary. The scales R+ and i?_ and the angle are 
determined by diagonalizing the potential Eq. ([13]). Because the short-distance wave functions ([15]) do not constrain the 
two components u{r) and w(r) as strongly as Eq. ([TT]) does in the case of OPE, any solution based on the Schrodinger 
equation ^ and the long-distance boundary conditions Eq. ^ can be matched to the r <C i?_ , i?+ wave function. 
Apart from an overall normalization, the long-distance wave function of the deuteron depends on two parameters, 
7 and rj. Once these parameters have been fixed, e.g. by experiment, we are able to calculate the wave function to 
arbitrarily small distances r > uniquely and regulator-independently. 

The numerical solution can thus be used to determine the parameters of the short-distance solution (I15p . In this 
way 7 and 77 are linked to two of the contact interactions that are present in the NNLO xPT potential. The third 
contact interaction does not affect the bound state, but can be determined by an examination of continuum wave 
functions. For example, for Set IV of the LECs we have i?+ = 2.11 fm and i?_ = 1.16 fm. Choosing the normalization 
of the long distance wave function such that u(r) —^ e~''^, and fixing 7 to its experimental value, the dependence of 
C-f and 0+,- on i] can be computed. The expressions can be found in Ref. [21j . where an explicit expression for the 
dependence of As on 77 for the case of this potential is also given. 

Since the dominant singularity of the NNLO xPT potential is r^^ the solution of the '^Si-'^Di bound-state problem 
when this potential is used for all r leads to a deuteron wave function that approaches zero more rapidly as r — > 
than does the LO xET wave function. As a consequence, deuteron matrix elements computed with wave functions 
where the NNLO potential is treated non-perturbatively have improved ultraviolet convergence properties |48| . In 
particular, matrix elements such as (1/r) and (l/r2), which are somewhat sensitive to short-distance pieces of the 
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r(fm) 



Figure 2. Comparison of radial deuteron wave functions u{r) and w{r) computed with the NNLO xET potential using a 
momentum-space regularization and the co-ordinate-space regularization of Ref. [21j . 



wave function, take values that agree quite well those found when sophisticated NN potential models are employed 
for their evaluation. 

Now we turn our attention to momentum-space solutions of Eq. ([2]) with the NNLO TPE potential derived in 
Refs. [28}i29j. We adjust the momentum-independent contact interaction Cs + Ct so that the experimental deuteron 
binding energy is reproduced. A numerical Fourier transform allows us to determine 77 as a function of the cutoff A. 
We find that for cutoffs A « 3.9 fm~^, 8.45 fm"'^, 14.1 fm~^ and 20.9 fm~^, the central value of the experimental 
D/S ratio, rj = 0.0256 [45,, is reproduced. Therefore at these values of the cutoff the coefficient of the second contact 
interaction needed to solve the ^Si-'^Di bound-state problem can be chosen to be zero. 

When the momentum-space problem is solved in this way and the momentum-space wave functions are (numeri- 
cally) transformed to co-ordinate-space the long-distance behavior matches that obtained in Ref. [21]. Therefore, for 
these cutoffs, we should find equivalent solutions in momentum and co-ordinate-space in the region where f ^ yx- 

In Fig. [21 we compare the Fourier-transformed wave functions of the momentum-space calculations and the 
regulator-independent result of the co-ordinate-space calculation. For the cutoff range that we consider here, we find 
agreement down to radii r « 1 fm. We hypothesize that \i A^ oo the wave functions agree for all r > 0. But as before, 
the momentum-space cutoffs need to be surprisingly large to find good agreement. E.g., for cutoff A k, 3.9 fm~^, there 
are visible deviations to the r-space wave function at r w 3 fm, and even for the largest cutoffs we observe deviations in 
the peak structure of w{r) around r « 1 fm. In the next section we will test whether this deviation has any noticeable 
effect on deuteron electromagnetic form factors. 



3 Deuteron form factors at leading order in ^ET 

Elastic electron-deuteron scattering can be parameterized by three independent form factors [I1[H[3], which we will 
here write as Breit-frame matrix elements of the two-nucleon four-current J^: 

Gc = ^ ((1 1 J"| 1) + (0 I J0|0) + (-1 1 J"| - 1)) ' (16) 
Gm = -^^(1 1 J+|0), (17) 



/2^|e 
1 

The form factors defined in Eqs. ([TB)) ([Hj) are related to the static moments of the nucleus by 



GQ = ^^i72((0k»-(lk1l))- (18) 



Gc(0) = 1, (19) 
GQ(0) = Qd, (20) 

Gm(0) = M.^, (21) 
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with Md the deuteron mass, Qd the deuteron quadrupole moment, and the deuteron magnetic moment in units 
of nuclear magnetons. In Eqs. (|16 |) -(fT8 |) . we have labeled the deuteron states by the projection of the deuteron spin 
along the direction of the three-vector p'g — Pe = q, and rj = (5^/(4Mj), with = |qp since we are in the Breit 
frame. We can then calculate the deuteron structure functions: 



A = Gl 



B 



(22) 
(23) 



In terms of A and B, the one-photon-exchange interaction yields a lab-frame differential cross section for unpolarized 
electron-deuteron scattering [2] 



da 



CTMott 



1 



^sin^ 



(^) 



(24) 



Here 9e is the electron scattering angle and E the electron energy, and CTMott is the Mott cross section. Measurement 
of the differential cross section only yields information on Gm and a combination of Gc and Gq. A third observable 
(usually t2o, a tensor polarization observable which is sensitive to Gq/Gc) must be measured if all three form factors 
are to be disentangled, and electron-deuteron scattering realize its full potential as a tool for measuring the deuteron's 
four-current. 

If we make a chiral expansion for the deuteron current operator J24','25','26J the results, up to corrections suppressed 

can be written as: 



by three powers of the xPT expansion parameter P = 



(pVo(q)lp) = |e|4^)(Q2)5(3)(p'_p_q/2), 



(p'|J(q)|p> 



|P 



q/4. 



G''^\Q-') + ifis<TxqGt/iQ 



M 



(5(3) (p' 



q/2), 



(25) 
(26) 

IVI 

(s) (s) 

with and GjJ the isoscalar form factors of the nucleon, and fis the isoscalar nucleon magnetic moment. In 
obtaining these equations we have dropped corrections to that have coefficients We will briefly discuss such 

"relativistic" corrections to the deuteron current operator in Sec. [4] 

Sandwiching the operators (l25|) and ([26|) between deuteron states yields the following co-ordinate-space integrals 
for the deuteron form factors [T]: 



2M 



Gc(Q') 
Gm{Q^) 



G^^\q') I dr[u\r)+w'{r)] JO , 



u{r)w{r) — 



J2 



\q\r 



(27) 
(28) 
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-32 



+G^)(g2)2 / dru\r)j, f 



+Gif(Q2) jV2 / dTu{r)w{r)j2 



qr 



— / dvw'^{r) 



Jo 



|q|r 
2 



J2 



|q|r 
2 



(29) 



where jo and j2 are spherical Bessel functions. 

We have demonstrated in Sec. [2] that we can determine regulator-independent deuteron wave functions in co- 
ordinate space. They may be used to compute the integrals Eqs. ((27|) -([29 |l . At the same time, for any finite value of A, 
we can evaluate the form factors Gc, Gq, and Gm using the wave functions obtained by solution of Eqs. ^ and ([2]). 

In the following, we present the results as a function of |q| = \/Q^ following the convention of Ref. 26J. We restrict 
ourselves to |q| below 1 GeV. This should cover essentially the whole range where xPT is expected to converge. 

We do this within the context of a chiral expansion for ratios of deuteron and nucleon form factors, e.g. Gc/G^\ 
so that we can focus on the predictions of xET for deuteron structure. Data for the single-nucleon form factors could 
then be employed to generate final results for deuteron form factors. But, up to the chiral order to which we work. 



this procedure is equivalent to using phcnomcnological input for G^' and G^^^ in the evaluation of expressions for 
Gc, Gm, and Gq, e.g. those in Eqs. ([27 |l ~([29 |) . Therefore we simply use the phenomenological single-nucleon form 
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Figure 3. Charge form factor of deuterium, Gc, obtained using different regularizations of the leading-order xEV potential. 
The red dotted line, green dashed, and blue dot-dashed lines are for momentum-space calculations at /l = 3, 10, and 20 
fm~^, respectively. The black solid line is the co-ordinate-space result. The experimental data is taken from the compilation of 
Ref. [51) : upward triangles represent data from the T20 measurement of Ref. [53], open circle [54], solid circle |55j . open squares 
[56] . downward triangles [57], rightward triangles [52], star [58], solid squares [59], solid diamonds [60] . 

factors of Ref. [49] in our evaluation of the deuteron form factors. These single-nucleon form factors are based on 
dispersion relations, and have the practical advantage of quantifying error bands due to uncertainties in the input 
experimental data. If we instead use a strict xPT expansion for the nucleon form factors our predictions for deuteron 
form factors begin to deviate from data for |q| ~ 300 MeV [25[|26j . However, this deviation has nothing to do with 
deuteron structure. Instead it is associated with the restricted range of |q| over which xPT provides a valid description 
of nucleon structure. This is the reason that we employ the single-nucleon form factors of Ref. [49j as an input in our 
calculation. Such a procedure is consistent as regards our xET examination of electron-deuteron scattering because 
deuteron and nucleon electromagnetic structure factorize in this reaction up to corrections that are O(P^) [26j . 

We first look at Gc based on the OPE wave functions: see Fig. [3] For A 00, the co-ordinate- and momentum- 
space wave functions yield essentially the same result. The approach to the yl — s- 00 limit is also shown in Fig. [3] 
and is quite interesting. At vl = 10 fm~^ the form factor has converged to the yl ^ 00 result. However, the result 
for yl = 3 fm^^ predicts a significantly different position of the zero of Gc- As we shall see below, this feature is 
particularly affected by higher-order corrections. The counterterm that would remove this yl dependence is one such 
effect, although it is of very high order: O(eP^). Small contributions of higher order will have most pronounced effects 
at such a zero, where all lower-order effects cancel by definition, and so sensitivity to the regulator scale in the vicinity 
of the minimum is not surprising. We note that for |q| < 600 MeV, the difference between yl = 3 fm~^ and yl = 20 fm~^ 
results is never more than a few per cent. 

The results obtained in co-ordinate-space and using large cutoffs in momentum space also agree very well for 
Gq. As seen in Fig. [4] convergence of the momentum-space results to the yl — > 00 limit is even more rapid there. 
Presumably, this is because, as was observed in Ref. [50], Gq is less sensitive to short-distance details than is Gc- There 
is, however, one exception to this: the value of Gq at = 0, i.e. the deuteron quadrupole moment, has significant 
evolution with yl, dropping by more than 10% between yl = 3 fm~^ and yl = 00. This can be seen in Table [21 where, 
for completeness, we also compile some other basic properties of our deuteron wave functions. As was observed in 
Ref. [27| . such yl-dependence is symptomatic that the O(eP^) counterterm impacts Qd more significantly than one 
would naively expect. A counterterm of natural size at 0{eP^) can absorb this large a cutoff dependence, and once it 
does, the remaining short-distance effects in Gq are minimal below |q| = 600 MeV [27] . 

In Gm the 00 limit again yields the co-ordinate-space result. But the results for low cutoffs deviate from the 
yl ^ 00 limit more than in any of the other form factors, as displayed in Fig. [5] They differ by 10% at |q| — 600 MeV, 
which corresponds to probing momenta of |q|/2 = 300 MeV in the nucleus itself. This sensitivity is entirely consistent 
with the presence of a counterterm for Gm at 0{eP^), which is only N^LO. Indeed, our calculation of Gm includes 
all NLO effects in the deuteron three-current J (apart from small 1/M^ corrections) so the counterterm would appear 
in J at the next order beyond what is presented here. 
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Figure 4. Quadrupole form factor of deuterium, Gq , obtained using different regularizations of the leading-order xET potential. 
Legend for theory curves as in Fig. [S] Experimental data taken from Ref. [51j . 



Table 2. Cutoff and scheme dependence of the kinetic energy T in MeV, the deuteron D-state probability Pd in %, the rms 
radius rd in fm, the quadrupole moment Qd in fm^, the asymptotic S-state normalization in fm~^''^ and the asymptotic S/D 
ratio 7]. The binding energy was fixed to -B = 2.225 MeV in all cases. For wave functions obtained in momentum space the 
cutoff A is given in fm"'^ while the r-space cut-off Vc is in fm. The first two sets of results are based on OPE, while the lower 
rows are results for TPE (with Set IV c'^s). Experimental values for observable quantities are listed in the last row, together 
with references. 
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Figure 5. Magnetic form factor of deuterium Gm obtained using different regularizations of the leading-order xET potential. 
Legend for theory curves as in Fig. [S] Experimental data from deuteron magnetic moment, open triangle [61]; the parameteri- 
zation of Ref. [3], open squares; and measurements of B(Q^): solid circles [62], open diamonds [63], and stars [64]. 



These issues aside, we see that the LO result for all three form factors is in remarkably good agreement with the 
data. In fact, as A —^ oo, the agreement with the data improves consistently. Obviously, the unphysical, spurious bound 
states affect the structure of the wave functions only at distances which are not relevant for momentum transfers |q| 
below 1 GeV. The spurious bound states in the '^Si-'^Di channel have no impact on these observables in the domain 
of validity of xET. Consequently their appearance in the spectrum in no way signals a breakdown of this approach. 

The agreement obtained between the data and the form factors evaluated with LO wave functions is surprising and 
also encouraging. One advantage of xET over phenomenological models of deuterium is its ability to systematically 
improve its predictions. Therefore in the next section we examine what happens when certain sub-leading effects are 
included in the wave functions and also discuss the impact on the predictions of pieces of the current operator that 
are of higher chiral order. 



4 Impact of two-pion exchange and exchange currents on the form factors 

To begin our discussion of higher-order effects on the deuterium form factors we first discuss some of the higher-order 
contributions to the current operators. Both Jq and J have a chiral expansion, and to produce the results displayed 
in the previous section we have ignored terms of order (^) . We have also dropped terms suppressed by (-p-) .In 
general corrections (^)" will be larger than corrections (ji)"^, as emphasized in Weinberg's original power counting 

used, e.g., in [28U66] . (^)" modifications to Eqs. ([27)) - ([29]) are associated with two-body operators because we have 
used phenomenological nucleon form factors. The strength of these terms is not fixed by any symmetry and must be 
determined by data (from either the single- nucleon or two-nucleon sector). While the first such effects in Gm are of 
relative order P^, the first occurrence of such a two-body operator in Gc and Gq is not until relative order P^. It is 
therefore higher-order than all the TPE effects in V that we will consider here. 

On the other hand, corrections to Eqs. (|27|) and ([28]l of nominal size [jj) occur in Jq with n = 2. These are 
corrections to the one-body part of Jq, and they can be obtained by demanding that matrix elements of the one-body 
current operator obey the correct transformation properties for Lorentz boosts to frames with velocities v <^ c. If we 

do this, we find a modification to the one-body part of Jq that introduces an additional factor of 1 — in Eq. ([25|l . 
In addition, if we are to have low-energy Lorentz covariance for deuteron matrix elements of the current operator 
(Jo, J) the momentum transfer |q| that is used in evaluating the integrals of Eq. ([27|l -(f29 |) must also be redefined. 
(This accounts for the boost of the deuteron wave function from the deuteron center-of-mass to the Brcit frame [SS] 
[26].) Both of these are small effects at the range of values considered here. Specifically, the former does not affect 
the minimum of \Gc \ at all, while the latter shifts it to the right by 3%, i.e. about 15 MeV for the LO results presented 
above. A two-nucleon effect in Jq that is of order tjtjj will be discussed explicitly below. 

Because these effects to the operator Jq are generically small, in what follows we focus on the impact of 

improving the wave functions. In particular, we first want to apply the NNLO wave functions that were discussed 
in Sec. [2] and which incorporate corrections of size (^) . The procedure used to evaluate Gc, Gq, and Gm, with 
these wave functions is the same as in Sec. [3] i.e. single-nucleon form factors of Belushkin, Hammer, and MeiBner are 
incorporated in the calculation to account for the effects of single-nucleon structure as per Eqs. ([27 |) -([29 | . 
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Figure 6. Quadrupole form factor of deuterium Gq obtained using using different regularizations of the TPE potential. Legend 
as in Fig. |3l with slight differences in cutoffs indicated. 
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Figure 7. Charge form factor of deuterium Gc obtained using using different regularizations of the TPE potential. Legend as 
in Fig. O with slight differences in cutoffs indicated. 



The pattern of convergence to the yl — > oo limit is very much the same as in the case of the OPE potential discussed 
in Sec.[3l Here, we first show results for Gc and Gq, in Figs. [6] and [71 respectively. Gq remains largely insensitive to 
A. At |q| = 0, the sensitivity is as large as that at any |q|, and is approximately 3% over the range from A — 3.9 fm^ 
to yl = 14.1 fm~^. Perhaps not coincidentally, this is roughly the size of the correction needed to move the value of 
Qd obtained with NNLO xET to the experimental result [27]. We find that Gc is more sensitive to A. In particular, 
the cutoff dependence of Gc around its zero is somewhat larger than naive expectations given the breakdown scale of 
xPT, and so we stress that that such expectations might be misleading, especially in the vicinity of a zero like this. 
For both form factors, the limit as yl — > oo of the results with the momentum-space cutoff is equal (within numerical 
uncertainties) to the result obtained via imposing a boundary condition in co-ordinate space at r = 0. 

In Fig. [Sjwe also show results for one of the modern xPT TV TV interactions of Ref. For this observable, we will 
show below that the uncertainties due to different sets of ttN LECs, the nucleon form factor and different choices for 77 
are rather small. But Fig. [6| shows that results based on the NNLO xET wave functions of Ref. [66], which employ a 
regulator exp{—p^/A^) with A in the range 500-700 MeV, do not overlap with the predictions from our wave functions 
at any |q|. 

Having demonstrated the equivalence of co-ordinate- and momentum-space formulations of TPE in the limit in 
which the regulator is removed, we now examine the impact of TPE on the form factors themselves by comparing 
results for the various sets of q's given in Table [TJ We begin with results for Gc, which are presented in Fig. [3 where 
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Figure 8. Comparison of charge form factor Gc with (blue, green, and orange broken curves) and without (black solid curve) 
the NNLO corrections included in the A'^A'' potential. The green double-dot-dashed, blue double-dash-dotted, and orange short- 
dashed curves use the different choices of the Ci's listed in Table [1] The error bands in the theoretical calculation shown for 
the "Set IV" choice of ttN LECs incorporate the experimental uncertainties in both the input value of 77 and the nucleon form 
factors of Ref. [49]. Data from the compilation of Ref. [51j . 



we also compare to data. Here we display the results for Sets I, III, and IV. (The results for Set II are very similar to 
those obtained with Set I.) The error band shown for the Set IV choice incorporates both uncertainties in the deuteron 
wave function (mainly those from varying the asymptotic D / S-raXio in this case), as well as the uncertainty bands in 
the single-nucleon form factors quoted in Ref. [49] . Sensitivity of Gc to these uncertainties is similar for the other two 
sets of ttN LECs. We note that such uncertainties are smaller than the sensitivity to different choices for the q's. 

The results clearly show that TPE corrections are small, as anticipated from the chiral expansion, and that — 
independent of the ttN LECs used to compute TPE effects — they shift the minimum of \Gc\ to the left. This yields 
a rather remarkable agreement between the Gc obtained once TPE is included in the deuteron wave function and 
experimental data. We will see below that effects of similar size can be expected from exchange-current contributions 
to the operator Jq. Therefore, at this point, no further conclusions on the size of LECs in TPE can be drawn. However, 
the results clearly indicate that it will be interesting to study Gc again taking such two-body effects in Jq into account. 
The results of Fig. [5] also suggest that the minimum of \Gc\ could be a fruitful place to look for improvements in the 
eft's description of data due to the inclusion of explicit /i(1232) degrees of freedom. Such calculations are, however, 
beyond the scope of this study. 

If the NNLO potential is used to compute deuteron wave functions the deuteron's asymptotic D-to-S ratio, ?/, is 
a free parameter of the calculation. In Fig. [5] we show how the position of | Gel's minimum depends on 77. Note that 
in order to make such a plot we need to specify the long-range potential, i.e. the set of q's we used for TPE. Also 
shown in Fig. [9] is the regulator-independent result obtained with the OPE potential, where 77 is not a free parameter, 
but is determined by the dynamics. This figure demonstrates that the position of the \Gc\ minimum is only weakly 
dependent on 77. In particular, had we used 77ope(7) ~ 0.026333 as our input instead of the central experimental value 
jyexpt = 0.0256 the minimum would only have shifted about 4% to the right. But the overall shift in the minimum 
from the OPE result to the TPE result with Set IV q's is 12% (the shifts for Set I-III q's are even larger). So roughly 
two-thirds of the shift in the \Gc\ minimum arises from the inclusion of TPE and corresponding changes of the wave 
function at distances r ^ Xjm-^. When regulator- independent wave functions are employed two-pion exchange has a 
significant impact on the |q| at which the zero of Gc occurs, regardless of which q's are used in its evaluation. 

Turning our attention now to Fig. [TOl we see that the OPE and TPE results for Gq are very close together. Above 
|q| = 350 MeV the effect of TPE is to shift Gq downwards, irrespective of the set of q's chosen. However, the error 
bands on the OPE and TPE calculations are always overlapping (recall that though we only show the error band 
for the TPE: Set IV calculation, all three TPE calculations have bands of similar size due to uncertainties in their 
input). The only place where an unambiguous difference between the different NN potentials can be seen is in the 
quadrupole moment. But here it is difficult to draw any definitive conclusion, since all calculations have central values 
that under-predict the experimental result = 0.2859(3) fm^ [671168]. This discrepancy is, however, consistent in all 
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Figure 9. Tlie position of tiie zero in the charge form factor, go, as a function of the asymptotic D-to-S ratio r] for the four 
different sets of ttN LECs considered here. The (Icr) experimental constraints on rj and go are indicated by the shaded region. 
The result found for r) and go from the LO xET potential with yl — + oo is given by the square. 
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Figure 10. Comparison of quadrupole form factor Gq with (blue, orange, and green broken curves) and without (black solid 
curve) the NNLO corrections included in the A'^A'^ potential. The green double-dot-dashed, blue double-dash-dotted, and orange 
short-dashed curves use the different choices of the d 's listed in Table [1] The error bands in the theoretical calculation shown 
for the "Set IV" choice of ttN LECs incorporate the experimental uncertainties in both the input value of rj and the nucleon 
form factors of Ref. |49) . Data from the compilation of Ref. [51) . 



cases with the expected short-distance contribution to Gq which occurs at 0{eP^) in the chiral expansion for Jq. All 
of this suggests that Gq is not a good place to attempt to test the TPE contribution to the deuteron wave function. 

In contrast, Gm is quite sensitive to the choice of dynamics for the long-range part of the NN potential, as is 
seen in Fig. 1111 The uncertainties from the input value of 77 and the single-nucleon form factors are sizeable, but even 
allowing for these uncertainties there are differences between OPE and TPE wave functions at the l~2a level. And at 
least at |q| > 0.5 GeV, there is also significant sensitivity to the choice of ttN LECs that is employed in the sub-leading 
TPE. Unfortunately the effect of TPE here is to worsen the excellent agreement with data that is achieved with OPE. 
(Although we note that Set IV is still marginally consistent with data if we allow for all input uncertainties.) However, 
it should be remembered that the theoretical uncertainty due to higher-order corrections in Gm is not depicted here. 
Such effects in the current J are suppressed by / relative to leading, the same relative size as TPE, and so 
should be included if a consistent calculation up to that order of Gm in the chiral expansion is desired. Until such a 
calculation is completed only tentative conclusions about the impact of the 0{P^) pieces of the NN potential on Gm 
can be drawn. 
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Figure 11. Comparison of magnetic form factor Gm with (blue, orange, and green broken curves) and without (black solid 
curves) the NNLO corrections included in the A'^A'^ potential. The green double-dot-dashed, blue double-dash-dotted, and orange 
short-dashed curves use the different choices of the Ci 's listed in Table [1] The error bands in the theoretical calculation shown 
for the "Set IV" choice of ttN LECs incorporate the experimental uncertainties in both the input value of t) and the nucleon 
form factors of Ref. |49) . Data as in Fig. [S] 

Similar caution is advisable in the interpretation of our results for Gc and Gq. In this case there is a two-body 
piece of Jo at relative order whose coefhcient is fixed by the low-energy consequences of Lorentz covariance. This 
operator was first derived in Ref. fTiT , and was re-derived in xET in Ref. [21? . The two-body matrix element of this 
piece of Jq can be computed reliably with the momentum-space wave functions of Section [2l It converges to a definite 
result as yl ^ oo and shifts the minimum of \Gc\ to the left by about 100 MeV, worsening the agreement with data 
seen in Fig. [S] However, as emphasized by Friar [71 and Adam and Arenhovel [72], this correction to Jq is associated 
with P'^ /A'P corrections to the OPE potential used in generating the wave functions. None of those P"^ /M"^ pieces of 
V were included in our calculations. Therefore a full evaluation of the relevant 1/M-suppressed effects remains to be 
performed with these wave functions. Such an evaluation is, however, beyond the scope of this paper. 



5 Conclusion 

Electron-deuteron scattering provides a window onto NN dynamics that gives complementary information to that 
obtained in NN scattering. It allows us to work at fixed total energy of the NN system and examine the electromagnetic 
response of the system as a function of the momentum-transfer squared, (which is equal to — |qp in the Breit frame). 
In xET the predictions for the three electromagnetic form factors Gc, Gq, and Gm that determine all elastic electron- 
deuteron scattering observables in the one-photon-exchange approximation are given solely by single-nucleon operators 
up to corrections of relative order P^ . Consequently the venerable formulae P71) - ([25)l for these form factors as Bessel 
transforms of probability densities are valid to this level of accuracy. 

Examination of these formulae indicates that the charge form factor is sensitive to distances up to a value of r 
given by 



where the factor of tt upstairs arises from the first zero of j'o(x) and the factor of 1/2 in the momentum transfer occurs 
because only half of q is transmitted to the relative degree of freedom. Thus as we move through the momentum 
range < |q| < 1 GeV Gc changes from an average of the probability density over all r to one where the important 
dynamics is that taking place at r w 1.2 fm. This radius is well inside that at which OPE is active, and indeed is 
small enough that some TPE effects are probed. Once |q| « 1 GeV we expect that two-body operators which scale 

as with n > 3 will no longer be suppressed, and so the dominance of the single-nucleon contributions to the 

deuteron four-current is no longer guaranteed there. But the momentum range |q| < 1 GeV already provides a wide 
kinematic domain over which we can test xET predictions for the deuteron wave functions (or equivalently — up to P^ 
corrections — form factors), and examine the impact of regularization on these predictions. 

Several calculations of these form factors that use xET already exist in the literature [24l[25l[26ll27j . But here we 
focused particularly on the possibility to obtain regulator-independent predictions for Gc, Gq, and Gm in the regime 
|q| < 1 GeV. We found that the results for all three form factors become stable as the momentum-space cutoff on the 
xET calculation, yl, is taken to infinity. The resulting predictions are therefore free of artifacts due to the particular 




(30) 
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function that is chosen as a short-distance regulator. They can also be obtained by employing a co-ordinate-space wave 
function where the Schrodinger equation is solved for all r > using NN potentials derived from chiral perturbation 
theory. We found that the momentum-space cutoff must be taken to be more than 2 GeV before stability is obtained 
in some cases (e.g. Gc near its zero, Gm at |q| « 0.45 GeV). We therefore caution that invoking naive-dimensional- 
analysis estimates to set the size of cutoffs that should be employed in xET can lead to the presence of significant 
cutoff artifacts in the results for observables. 

We were able to perform these calculations for both the LO NN potential (which consists solely of one-pion 
exchange at long range) and for the NNLO NN potential. Because both of these potentials are singular and attractive 
as r ^ they can generate a shallow (deuteron) bound state. But their singular nature means that as A oo they 
also generate spurious bound states at energies beyond the range of applicability of the theory. However, these spurious 
bound states have no impact on the deuteron form factors, and we found no indication that their appearance signals 
a breakdown of the xET approach. 

We were not able to perform a calculation with the NLO xPT NN potential, since the potential at this order is 
singular and repulsive. In consequence it is impossible to take the yl — > oo limit with this V inserted in the Schrodinger 
equation and still have a spectrum with a shallow bound state. An order-by-order assessment of the convergence of 
xET must therefore await the development of a power counting which does not have this deficiency. In this context 
we note that the version of xPT in which the Z\(1232) is included as an explicit degree of freedom leads to a potential 
at NLO in the (modified) chiral expansion that is singular and attractive. Indeed, in the '^Si-^Di channel the NLO 
potential in the theory with explicit Deltas is very similar to the NNLO potential we have used here. Therefore we 
believe that a study employing the NLO V from xPT with explicit Deltas — while definitely called for — will likely 
reach similar conclusions to those found here using the NNLO V from standard (Delta-less) xPT. 

The calculations to LO in xET gave predictions for deuteron electromagnetic structure which are remarkably close 
to experiment for |q| < 0.6 GeV once the yl ^ oo limit was taken. The only significant disagreement between the LO 
prediction and the data in the compilation of Ref. [FT] occurred in the vicinity of the zero of Gc- But even there the 



The corrections to the LO result once the NNLO potential was employed to compute deuteron structure were also 
consistent with this expectation, proving to be small, especially in the case of Gq, which is already well described by 
the LO calculation. Our results also show that the discrepancy between the LO xET calculation and experimental data 
for Gc in the vicinity of its zero may be eliminated when such a higher-order calculation is performed. In particular, 
we showed that the position of the minimum in Gc is sensitive to details of the two-pion-exchange interaction that 
is present in the NNLO xPT potential. This will be an interesting point to focus on in future investigations of two- 
pion exchange. In contrast, the good description of Gm in LO appears to be accidental. Various small but significant 
higher-order effects apparently cancel for this observable, but this needs to be checked in a complete higher-order 
calculation. 

We also assessed the impact of the 0{P^) corrections to Jq on our results for Gc and found that they converge to 
a definite limit as yl — *■ oo. Although we did not include any xPT corrections to the leading-order in the results we 
presented, our estimates of the 0{P^) effects in Jq show that the exchange-current corrections have roughly the same 
impact on Gc in the vicinity of its zero as do the 0{P^) contributions to V . This suggests that the power counting 
for the potential and the currents is working well in this domain. 

All of this makes it very interesting to complete consistent higher-order calculations of deuteron form factors using 
these regulator-independent wave functions. To do this the O(P^) pieces of will have to be computed fully and in a 
manner that is consistent with the treatment of corrections in V. Only when such a computation is performed 

will we know if the improvement found here with the NNLO X-ET deuteron wave functions represents genuinely good 
convergence of the chiral expansion for these observables or not. The success we observed here may be merely a 
fortunate result of examining only the corrections to V . And an assessment of the impact that 0{P'^) (and higher) 
pieces of Jo will also be necessary if a definitive conclusion as regards the impact of TPE contributions on the position 
of Gc's zero is to be reached. Irrespective of these issues, though, it is clear that the deuteron charge form factor is 
sensitive to details of the chiral dynamics that is at work in deuterium. 
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discrepancy between theory and data was within the expected 




, which is about 25% at 1 GeV. 
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